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Abstract 


The study of inclusive semileptonic decays of B mesons is analyzed from 
the viewpoint of probing hadron structure and strong interactions. General 
formulas for the differential decay rates are given in terms of the structure 
functions in arbitrary frame of reference, taking into account the finite charged 
lepton mass. These formulas can be used for structure function measurements. 
The behavior of the structure functions is shown to contain information on 
the constitution of the B meson and the dynamics of strong interactions. 
Measurements of the structure functions would reveal the existence of a point¬ 
like b quark in the B meson, establish the spin 1/2 nature of the b quark, and 
identify the b quantum number of the B meson. These measurements would 
determine the b quark distribution function, and allow new insight into the 
nature of confinement. The b quark distribution function can also be extracted 
directly through the measurements of the decay distributions with respect to 
the scaling variable 
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I. INTRODUCTION 


Coloured quarks and gluons as physical constituents of hadrons and quantum chro¬ 
modynamics (QCD) as the theory of strong interactions have been generally accepted as 
the modern point of view with regard to the structure of hadrons and strong interactions. 
Asymptotic freedom of QCD explains Bjorken scaling [Q] observed in deep inelastic electron- 
proton scattering before the advent of QCD and allows perturbative calculations of strong 
interactions at high energy or short distance. On the other hand, the increase of the strong 
coupling constant of QCD at low energy or long distance provides a reason to accept quarks 
and gluons as physical constituents of hadrons without directly detecting free quarks and 
gluons. There are a variety of experimental results, which not only provide conhrmations of 
the quark and gluon structure of hadrons, but also lead to incisive precision tests of QCD. 
While the evidence is compelling that QCD is a realistic theory of strong interactions, there 
is still a very incomplete understanding of the nonperturbative, conhning, long-distance sec¬ 
tor of QCD, which is a basic challenge in describing strong interactions, in particular how 
quarks and gluons bind together to form hadrons. 

The B meson is heavy and relatively long-lived . This makes it advantageous to probe 
hadronic physics. One powerful way of experimentally investigating the structure of the B 
meson and strong interactions is to study its inclusive semileptonic decay. The structure 
of the B meson and the strong interaction can be studied with this underlying weak decay 
process because the decay rate is modified by the strong interaction, which is responsible 
for quark conhnement and gluon radiation. The power and beauty of inclusive semileptonic 
decay is that the electroweak held generated during the semileptonic decay is well understood 
and leptons do not interact strongly. Furthermore, because of inclusive characteristics of 
the decay we would expect that the decay rate is mainly affected by the initial bound state 
effect and is insensitive to hadronization in the hnal states, which are summed over and 
possess the completeness. These permit us to probe the structure of the B meson and 
strong interactions by means of a known charged weak current. 

While much attention has been focussed on the determination of \ Vub\ and \ Vcb\ by inclu¬ 
sive semileptonic decays of B mesons, a wealth of information on the hadron structure and 
strong interactions available from these processes deserves intense scrutiny in experiment 
and theory as well. The purpose of this paper is to discuss the study of hadron structure 
and strong interactions by means of inclusive semileptonic B decays. We would like to show 
that this study will provide evidence for the quark structure of the B meson and allow new 
insight into the nature of conhnement. Of course, this study will also be benehcial to precise 
determination of the fundamental standard model parameters |Wb| and \ Vcb\- 

The outline of this paper is as follows. We derive in Sec. II within the framework of the 
standard electroweak model general formulas for the differential decay rates for inclusive 
semileptonic B meson decays in terms of the structure functions. The derivation is made 
in arbitrary frame of reference and the hnite charged lepton mass is taken into account. 
These formulas are useful for measurements of the structure functions. The behavior of the 
structure functions is explored in Sec. Ill, where we also elucidate by measurements of the 
structure functions what can be learned on the constitution of the B meson. The extraction 
of the b quark distribution function, to which the structure functions are related due to the 
light-cone dominance, is discussed in Sec. IV. In Sec. V we conclude. 
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II. STRUCTURE FUNCTIONS 


We want to compute the decay rate for the inclusive semileptonic decay of the B meson 
to lowest order in the weak interaction in arbitrary frame of reference. In this inclusive 
process B Xivi-, the B meson of mass M with four-momentum P decays into any 
possible hadronic hnal state X, over which we sum, via the emission of a virtual W boson 
which materializes as a charged lepton ^ of mass Mi with four-momentum ki and a massless 
antineutrino vi with four-momentum ky. By standard steps we hnd the unpolarized decay 
rate to be proportional to the leptonic tensor times the hadronic tensor W^y-. 


{2TTfE 2Ei 2E^ ’ 


( 2 . 1 ) 


where E, Ei, and E^ denote the energies of the B meson, the charged lepton, and the 
antineutrino, respectively. Vgb is the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, 
which is Vcb for the b ^ c transition induced decay and Wfe for the b —>■ u transition induced 
decay. The leptonic tensor for the lepton pair is completely determined by the standard 
electroweak theory since leptons do not have strong interactions: 


L"- = 2(k^k: + - J'-A,. A„ + 'if'*" „;,A,“A«). 


( 2 . 2 ) 


The hadronic tensor incorporates all nonperturbative QCD physics for the inclusive semilep¬ 
tonic B decay. It is summed over all hadronic hnal states and can be expressed in terms of 
a current commutator taken between the B meson states: 
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(2.3) 


where q = ki + k^ stands for the momentum transfer from the B meson to the lepton pair 
and j^iy) = g( 2 /) 7 ^(l — 'y 5 )b{y) is the charged weak current. The B meson state \B) satishes 
the standard covariant normalization {B\B) = 2E{27i)^6^{0). 

The most general hadronic tensor form that can be contructed is a linear combina¬ 
tion of PiiPu, P^iQu, qiiPv) diidu, ^iLuapP'^q^, and with coefficients being scalar functions 
W4(n, q^) of the two independent Lorentz invariants, v = q ■ P/M and q^. However, the 
combination P^qu — q/iPu does not contribute since {P^q^ — q^Py) = 0. Thus the hadronic 
tensor must take the form 


TIT TIT ^ ^ TIT TIT 

Wfyy = -g^yWi -F - ^£,dyag^^W3 + + 




M2 


-w 


(2.4) 


Equation (p.3|) shows that = Wy^j,, so Wa, a = 1,..., 5, are real. The interesting 
physics describing the hadron structure and strong interactions is wrapped up in the hve 
dimensionless real structure functions ITa(i/, q‘^),a = 1,..., 5, for the unpolarized processes. 

The unpolarized differential decay rate can be expressed in terms of the structure func¬ 
tions: 


d^T 

dgdudq'^ 


^gll^HV, 2(92 - M|) + - irf - + M}) 

+W 32 {uq‘^ - 2r]q^ + M/v)/M + W^Mf{q^ - M/)/M^ 
+WbAM/{v-g)/M], 


(2.5) 
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where we define the Lorentz invariant t] = ke ■ P/M. Equation (|2.5|) shows that the B 
meson structure functions could be separated and extracted from the measurement of the 
differential decay rate at different values of rj for fixed i> and q^. 

It would be useful to express the differential decay rate in terms of the angle 9 between 
the charged lepton and the B meson in the virtual W rest frame. The angle 6 is related to 
the Lorentz invariants p, u, and by 

/U'2 1^2 - 

2r, = {l + -^)u - (1 - -^)Ju^-q2cos9 . (2.6) 

qz qz V 

Using Eq. (PT 6 |), the differential decay rate in terms of cos 6 ', z/, and can be obtained 
from Eq. ( ^7^ : 

[W,2(q^ - Mf) + W,{v^ -f + Mj - 

+IU32(g2 - M2)^z/2-g2 cosO/M + W^M/{q^ - M/)/M^ 

+W52M/{u- n)/M], (2.7) 


d^T 

dcosOdvdq^ 


where 


Ml 


Ml 


K = —— (1- 2 “) V ~ 


( 2 , 8 ) 


From Eq. (^77|) we see that the structure functions could also be separated and extracted 
experimentally by looking at the angular 9 distribution of the charged lepton at each value 
of z/ and g^. 

Integrating Eq. {Tl) over cos 6 * gives 


Gl\V,i\\ Ml 

dvdq^ 327 i^E g2 ' 

+W^\\2vHi + ^- 

3 g2 


z/2 - g2{IUi2(g 
2Mf 


M|) 


g^ 


0 Mf 

qz 


Mf 


Mil t 2 T^/r2\ , TTA 2M|z/ M, 

- Mf)+^ 


J} 


(2,9) 


This decay distribution would also be useful to measure the structure functions except for 
IT 3 , whose contribution is absent. 

Equations (^.Sj), ( |2.7|) , and (|2.9|) show that the contributions of the structure functions 
IT 4 and lUs are suppressed as they are multiplied by the square of the charged lepton mass. 
Therefore, the inclusive semileptonic B decay with a tan lepton in the final state would be 
appropriate for measuring IT 4 and IT 5 . 

In the inclusive semileptonic B decay where an electron (or muon) is the final state 
lepton, the mass effect of it may be negligible. Neglecting the charged lepton mass, then the 
differential decay rates in Eqs. (|3p, (O), and (p.9|) reduce, respectively, to 


d^V 

drjdudq'^ 


32n^E 


[Wi2q^ + - W 


q^) + fha 


2g2 

~M 




2h)], 


( 2 . 10 ) 
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d^T 

dcosOdvdq^ 


GlW 


Y^z /2 — q^\Wi2q^ + — q'^)sin^ 6 * + W 3 cos^] , 


( 2 . 11 ) 


= + (2-12) 

Note that only the structure functions hhi and hh 2 contribute to the double differential decay 
rate d?T / {dudq^) with lepton masses being neglected. 

The formulas for the differential decay rates derived in this section are quite general and 
can be used to measure the B meson structure functions in any particular frame. However, 
we should point out that, in addition to the lowest order contribution in the weak interaction, 
there are higher order radiative processes which contribute to the decay rate. In order to 
extract the structure functions from experimental data, the radiative corrections to the decay 
rate must be implemented. The electroweak radiative corrections were discussed in Ref. [Q. 
The QCD radiative corrections were calculated at the parton level of quarks and gluons in 
Refs. 1^. 


III. BEHAVIOR OF THE STRUCTURE FUNCTIONS 


The behavior of the structure functions manifests the structure of the B meson and the 
dynamics of strong interactions. In this section we explore theoretically how the B meson 
structure functions look like. We will also discuss the previous results in Ref. 0. 

In order to see the way in which the structure of the B meson and strong interactions 
manifest themselves in the structure functions, it will be useful to start with an investigation 
of the behavior of the structure functions in the free quark limit when no strong interaction 
takes place. The unpolarized tree-level decay rate for the free point-like quark decay b qiu 
is given in arbitrary frame by 


GlWr u.d"ked^K 
(27r)'5pO 2E, 2E, ’ 


(3.1) 


where stands for the four-momentum of the b quark of mass mi,. is the leptonic tensor 
given in Eq. (|2.2|) . is the quark tensor 

= 4(5[(pb - qf - -q-Pb) + 2p^Pb + ie^"''^'"Pbpqa - {.pU'" + q^Pb)) > (3-2) 

where is the mass of the decay produced g-quark. 

In the free quark limit, the B meson and the b quark in it have the same velocity: 
Pb/mb = P/M. By use of the delta function relation 


where 
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z/± 




M 


(3.4) 


comparing Eqs. 
functions in the 


( p.l|) and (|2.4|) with Eqs. (p.l|) and ( p.2|) 
free quark limit 


allows us to identify the structure 


lyfree 

Wiree 

W^ee 

Wfree 




8^+ 


e+-e- 

4 


e+-e- 


5(e+- 


0 , 


mb . 

m’'‘ 

_ !!]± 
M 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 


The “free” functions g^)/2, (,^_,_ — .^_)hE|''®®(i/, g^)/ 8 , (.^+ — .^_)lT|( 5 ®(z/, g^)/4 exhibit 

the intriguing property that they are functions of only one variable and not of u and 
independently. It is interesting to see if this scaling property persists under the strong 
interaction modihcation, which we shall discuss below. 

We go on to explore the actual behavior of the structure functions as functions of their 
arguments u and switching on strong interactions. In the inclusive semileptonic B decays, 
the square of the momentum transfer lies in the range Mf < g^ < (M — where 

is the minimal hadronic invariant mass in the hnal state, which is the D meson (pion) 
mass for the b ^ c {b ^ u) transition. Because of heaviness of the B meson, the decays occur 
mostly at large momentum transfer. From Eq. ( p.3|) we see that at large momentum transfer 
the hadronic tensor is dominated by the space-time region near the light cone ^ 0 , where 
y is the space-time interval between the points at which the currents and ji/( 0 ) act. 

The light-cone dominance implies that the hve structure functions, a priori independent, 
are related to a single distribution function |]9HTT1: 


Wi 

W 2 

IE 4 

Wb 


2[/K+) + /K-)l, 

t ^, [(+/((+) - (-/((-)] > 

s+ ~ s- 

0, 


(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 


where the distribution function is dehned by l^-H 


m = / d(!/ ■ - 75)%)|B)I,».o . 


(3.15) 


The distribution function reflects the long-distance strong interaction responsible for the 
5-quark conhnement in the B meson. The structure functions show a simplihed structure at 
large momentum transfer. We see that on the light-cone dominance IT 4 remains zero and 
also IT 5 = IT 3 just as in the free quark limit. 
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It is instructive to pin down the form of the b quark distribution function in the free 
quark limit for comparison with its form in the real world. In this limit, the free Dirac held 
b{y) = with the velocity equality Pb/^b = P/M, so from Eq. ( |3.15|) it follows that 

the distribution function becomes 


/(O = 


mb . 

M' 


(3.16) 


Substituting Eq. ( b.l 6 |) into Eqs. (|3.10|) - (|3.14|) , we recover the structure functions in the free 
quark limit in Eqs. (|3.5|) - (|3.9|) . Deviations from the delta function behavior would signal 
the strong interaction effects. 

Several important properties of the distribution function can be derived from held theory 
PJT^. It obeys positivity with a support 0 < .^ < 1. The distribution function is normalized 


to unity: 


CdU{i) = i. 

JO 


(3.17) 


The normalization of f{^) is exact and does not get renormalized as a consequence of b 
quantum number conservation. Since the b quark inside the B meson behaves as almost 
free because of its large mass, relative to which its binding to the light constituents is 
weak, we expect that the distribution function is very close to its asymptotic form — the 
delta function in Eq. (|3.16|) . This is supported by the analysis using the operator product 
expansion and the heavy quark ehective theory method [|I3HT7||, which indicates that the 
distribution function is sharply peaked around ^ = rUb/M with a width of order Aqcd/M 

B. 

It is convenient to redehne the structure functions as 


10.12 


Fi{i^,q^) 


(3.18) 

F2{i^,q^) 


(3.19) 

F3{j^, q^) 


(3.20) 

F,iiy,q^) 


(3.21) 


With these dehnitions, it is straightforward to get from 


Eqs. (pTOD-fCT): 


Pi — F3 _ F2 + ^-Fi _ F2 + 

2 “ + ^ “ e+-e- 


/(e+)• 


(3.22) 


Note that on the light-cone dominance T 5 = T 3 . Equation ( ^. 22 |) shows scaling of the 
linear combinations of the structure functions, i.e., for large they depend only on the 
dimensionless scaling variable .^+. The normalization condition (|3.17|) of the distribution 
function implies the sum rules 



Fi — F 3 


2 



F2 + i-Fi 

i+ + i- 



F2 + ^-Fj 


1 . 


(3.23) 
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We may go further. Assuming the /(^-) terms in Eqs. (|3.1CI|) - (|3.14|) to be negligible 
implies the following relations of the structure functions |§] 


e+Fi = f 2 = -e+Fs = -e+Fs = e+/(e+) , 


(3.24) 


where the structure functions also scale, displaying the same property of the structure func¬ 
tions as in the free quark limit shown above. Furthermore, the following sum rules can be 
obtained from the normalization condition (|3.17|) : 


f mq + M^ 
M 


d^+F,= 


L+m, ^^^+^ 2 /.^+ = 

M 


mq+M£ 

M 


F, = 


-Ps = 1 ■ (3'28) 

M 


The relations ( p.22 ) and ( 3.24 ) are a consequence of the spin 1/2 nature of the b quark; 
spin 0 partons would lead to Fi = 0 or Wi = 0. 

The /(^-) term is a consequence of held theory, corresponding to the creation of quark- 
antiquark pairs in the hnal state. Here we pursue the arguments of 0 and examine more 
closely the assumption that the f{^-) terms in Eqs. (|3.1CI|) - (|3.14|) are negligible. The kine¬ 
matic ranges of and are 


rUg + Me 


M 


<e+<i, 


Me-mg^Me + < ^-< 1 


M 


M 


M ’ 


(3.26) 

(3.27) 


respectively, where Q{x) is a step function. For b ^ c decays ^ 0.5, the /(^-) terms 
entering Eqs. (|3.10|) - (p.l4|) can be safely neglected since the distribution function is sharply 
peaked near one. However, for 6 —> u decays can be as large as one, so that the sizes 
of /(^-) and /(^+) may be comparable. Thus the approximate scaling and the sum rules 
in Eqs. ( p.24| ) and (|3.25|) would hold for the inclusive charmed semileptonic B decay, but 
would be considerably violated by the non-negligible f{^-) term for the inclusive charmless 
semileptonic B decay. By contrast, the approximate scaling and the sum rules in Eqs. (p.22|) 
and ( |3.23| ) are valid regardless of the size of /(^-). 

Also worthy of noting is that there are other sources of scaling violation: logarithmic 
corrections due to gluon radiation and corrections of inverse powers of due to departures 
from the light cone. 


IV. EXTRACTION OF THE DISTRIBUTION FUNCTION 


The b quark distribution function f{^) dehned in Eq. (|3.15|) contains important infor¬ 
mation on the nature of conhnement. As a necessary input, its knowledge is also crucial 
for improving the theoretical accuracies on the determinations of \Vub\ and |V/b| from the 
charged lepton energy spectra |l0| , pd]| and the hadronic invariant mass spectrum |^. An 
experimental extraction of the distribution function will therefore be of special importance 
and interest. In this section we shall discuss how one could extract the distribution function 
from experiment. 

Equations (|3.10|) - (|3.14|) imply that /(^) can be extracted through the measurements of 
the structure functions, which is discussed in section H. Recently, it has been put forward 























































that a measurement of the differential decay rate as a function of can be also used 
to extract directly the b quark distribution function. The expression for dT/d^^ has been 
derived in []^, neglecting lepton masses. Here we generalize this result to the case where 
the charged lepton mass cannot be neglected. Substituting Eqs. (p.l0|) -( p.l4|) into Eq. ( p.9D , 
changing the variables from (z/, to and integrating over with the kinematic 

limits Mj < < {^+M — mq^, we arrive at 


dr 


Gl\Vqh\^M^ 

1927r3 E 


Cf{^+mx,z) + 


12 riUM-nig)^ 

Jm? 


dq^^v'^ - gV(^-)(l - 

V qz 

(4.1) 


with 


4 4 44\1 1+a^^ — + 


<h(a;, z) = y/X[l - 7(x^ + + / + x'^z‘^ + + 12x‘^z‘^ + x® + /] 

+12(2x"^ — z"^ — x"‘;2"‘)ln- 

+ 12(1 - x^); 2 ^ 1 n- 


2x 


4 4 {1 - X^f - {1 + X'^)Z^ + {1 - x‘^)^/\ 


2xz‘^ 


(4.2) 


n{i^,q^)=q^-Mj + 


2Mjv 

1^ 


, Mj, !£-, o o 


2Mf 


o, Mj 
4z/2 1 + -L 

qz 


where 


\ /l 2 |/i| 2 2\2 

A = -4x +(l + x -.) . ^ = - = 


2 + m 

(4.3) 


(4.4) 


z/ and are related to .^+ and by 


z/ = 


e- = 


^2 m 2 + g 2 - m 2 


2 e+M 

2 


m 


e+M 2 


(4.5) 

(4.6) 


If the charged lepton mass is neglected, Eq. ([4.1|) reduces to the expression given in . 


The contribution of the /(^-) term to the specihc spectra dT /d^+ in Eq. (|4.1|) is expected 
to be negligible since (1) the distribution function is very close to its asymptotic form, 
— rrih/M), in the free quark limit as discussed in the last section and in this limit f{^-) 
vanishes and (2) the dominant contribution to the f{^-) integral at a given results from 
the large region, corresponding to the neighbourhood of the upper integration limit for 
g2, which is suppressed by z/2 — g2. The numerical studies carried out in [|^ have conhrmed 


this expectation. The smallness of the /(^-) term has already been found in the parton 
model |2^. As a result, dr/d^+ is proportional to /(^+), so the measurement of the decay 
distribution can be used to extract directly the nonperturbative distribution function. In 
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other words, the decay distribution with respect to the scaling variable probes the QCD 
dynamics in a most sensitive manner. 

In practice, as the measurements of the structure functions, radiative corrections must 
be implemented in order to extract the distribution function from the measured raw yields 
using the above expressions. 

The decay distribution as a function of is a measurable quantity with many other 
remarkable uses. New hadronic uncertainty insensitive methods have been proposed 
for precise determinations of |I4b| and \Vcb\ through the measurements of the b u and 
b —> c decay distributions as a function of Large theoretical uncertainties associated 
with nonperturbative strong interactions can be avoided with these methods. Moreover, in 
these measurements rare 6 —>■ u events can be separated efficiently from 6 —>■ c events. These 
methods would enable us to determine especially \Vub\ in a model independent way. The 
b ^ u spectrum dT{B —>• Xuii')/d^+ could be also used to measure reliably the inclusive 
charmless semileptonic branching fraction of the B meson. 


V. CONCLUSIONS 

Inclusive semilpetonic B decay is a very powerful means of probing hadron structure 
and strong interactions. There are clear similarities between inclusive semileptonic B decay 
and deep inelastic lepton-nucleon scattering. Inclusive semileptonic B decay can be used to 
pursue the QCD as deep inelastic scattering. We have presented general formulas for the 
differential decay rates for the inclusive semileptonic decays of the B meson in terms of the 
structure functions with the inclusion of charged lepton masses. These formulas are suitable 
for measuring the B meson structure functions in any particular frame. 

Measurements of the B meson structure functions would provide evidence for the quark 
structure of the B meson. The scaling behavior of the structure functions at large momentum 
transfer as in Eqs. (|3.22|) and (|3.24|) would reveal the existence of a point-like b quark in the 
B meson. The experimental verification of nonvanishing Wi would indicate that the b quark 
has spin 1/2. The sum rules as in Eqs. (|3.23|) and (|3.25|) would provide a basic element in 
identifying the b quantum number of the B meson. 

The light-cone dominance simplifies the nonperturbative description of inclusive semilep¬ 
tonic B decays and permits the structure functions to be related to a single distribution 
function, leading to the testable theoretical predictions. The distribution function f{^) is 
defined as the Fourier transform of the matrix element of the bilocal and on-light-cone op¬ 
erator in 6-quark fields taken between the B meson states, which gives the probability of 
finding a b quark with momentum inside the B meson. This function depends only on 
the initial B meson state and is therefore universal. The nonperturbative distribution func¬ 
tion can be extracted by the measurements of either the structure functions or the decay 
distribution as a function of 

In addition to testing the theoretical predictions discussed here, advancing our under¬ 
standing of hadron structure and strong interactions, and providing new observations on the 
nature of confinement, measurements of the B meson structure functions and the universal b 
quark distribution function would be beneficial to precise determinations of the fundamental 
standard model parameters |I4fe| and |V/b|- 


10 












The study of hadron structure and strong interactions by exploiting inclusive semileptonic 
B decays will considerably extend the B physics reach. Measurements of the B meson 
structure functions and the b quark distribution function should prove rewarding. 
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